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ABSTRACT 

We investigate the prospects of future supernova searches to get meaningful constraints about the 
cosmic star formation history (CSFH) and the delay time of type la supernovae from star formation 
(ria), based only on supernova data. Here we parameterize the CSFH by two parameters, a and (3 that 
are the evolutionary indices [oc (l+z)^°^'^^] at z > 1 and < 1, respectively, and quantitatively examined 
how well the three parameters (a, /3, and Tia) can be constrained in ongoing and future supernova 
surveys. We found that the type classification of detected supernovae down to the magnitude of 
Iab ~ 27 is essential, to get useful constraint on (3. The parameter rja can also be constrained within 
an accuracy of ~ 1-2 Gyr, without knowing a that is somewhat degenerate with ria. This might 
be potentially achieved by ground-based surveys but depending on the still highly uncertain type- 
classification by imaging data. More reliable classification will be achieved by the SNAP mission. The 
supernova counts at a magnitude level of Iab or Kab ^ 30 will allow us to break degeneracies between 
a and ria and independently constrain all the three parameters, even without knowing supernova types. 
This can be achieved by the SNAP and JWST missions, having different strength of larger statistics 
and reach to higher redshifts, respectively. The dependence of observable quantities on survey time 
intervals is also quantitatively calculated and discussed. 
Subject headings: cosmology: observations — supernova:general 



1. INTRODUCTION 

It is important to measure the cosmic supernova (SN) 
rate at present as well as at high redshift, to understand 
the cosmic star formation history (CSFH) and to get a 
hint for the origin of the type la supernovae (SNe la). 
Although the CSFH up to z ~ 5 has been intensively 
studied based on star formation rate (SFR) indicators of 
high-z galaxies (Lilly et al. 1996; Madau et al. 1996; 
see Totani & Takeuchi 2002 for a recent compilation of 
observed data by various methods), the derived CSFH is 
still a matter of debate since the estimated SFR densi- 
ties are considerably different depending on authors and 
methods. The CSFH inferred from SN rate evolution in- 
cludes completely independent and complementary infor- 
mation. CSFH estimated by galactic SFR, especially at 
higher redshifts, is inevitably biased toward bright galax- 
ies; in many studies SFR was corrected by extrapolating 
the galaxy luminosity function to the highly uncertain 
faint end, well below the observational limit. However, 
SNe in any galaxy including dwarfs, or even those in truly 
intergalactic space, can equally be detected by SN sur- 
veys and hence it is completely free from this bias. The 
galactic SFR indicators and core-collapse SN rate gener- 
ally trace the massive star formation, but the SN la rate 
gives some information for low or intermediate mass (typ- 
ically a few solar mass) star formation (e.g., Yungelson & 
Livio 2000). The delay time of SNe la from star forma- 
tion to the explosion may be inferred from observations, 
giving a useful hint for the still unknown progenitors of 
SNe la. 

There are a number of publications about theoretical 
predictions of high-z S N rate as a probe of CSFH and the 
progenitor of SNe la ijTotani et alJ Il996t iMadau et all 
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these papers, only the time evolution of the cosmic SN 
rate density, which is simply converted from CSFH, was 
shown. However, to compare with observed data it is 
necessary to predict the expected number of SNe lim- 
ited by a given sensitivity of an obs ervation. A work in 
this direction was first presented bv lDahlen fc FranssonI 
(1999, hereafter DF99). 

On the other hand, in recent years there was 
a remarkable progress about measurement of type 
la and core-collapse supernova rate, pushing out 
the maximum distance from z ~ t o ^ 1.0 
(Painctal. 1996; CaDoellaro et al. 1999; Har din et^J 
(2000; Gal- Yam ct al. 2002; Madgwick et al. 2003 
iGal-Yam k. Maoz. 2QM. Dahlen et al.. .20041 . Some 
interesting implications on the CSFH and the la 
delay time have already been obtaine d by analyses 
of these data (Pain ct al. 2002; Gal- Yam et alJ 1200^ 
Maoz fc G al- Yam 2004; Blanc et al.ll2004HStrolger et alJ 



20041) ■ though statistics is not large enough to derive 



decisive conclusions. For example, one must use infor- 
mation on CSFH derived by other methods than SNe to 
derive a strong constraint on the la delay time. Future 
SN surveys will provide much better constraints and 
may allow to determine the interesting parameters only 
by supernova information. 

However, the cosmic SN rate evolution is dependent on 
a number of parameters that are often degenerate with 
each other, and it is still unclear how to get useful in- 
formation from future SN surveys. The purpose of this 
paper is to examine what can be learn from future SN 
surveys, and find the best strategy of analyzing data to 
get useful information. We set our goal as to constrain 
the three interesting parameters only based on informa- 
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tion obtained by supernova surveys: the rate of increase 
of the cosmic SFR from z = to ~ 1, the evolution of 
the SFR beyond z ~ 1, and the SN la delay time, by 
breaking degeneracy between these parameters. We will 
make predictions on various observational quantities for 
future surveys and discuss how well these parameters can 
be constrained in the near future. 

SNe are detected by variability in more than two 
epochs of observations. If a reference frame that is tem- 
porally well separated from new observations is available, 
we can measure the real SN flux relative to the refer- 
ence frame. However, the first SN candidates are often 
selected from a flux variability on a shorter time base- 
line that is comparable to the expected SN time scale 
(i.e., month), to find SNe at early brightening phase, and 
to avoid contamination of active galactic nuclei (AGNs) 
having longer variability time scales. The latest deep SN 
survey by GOODS was also based on a short time vari- 
ability. This selection process should have some effect on 
the event rate and hence on SN rate studies, even though 
the real flux and light curve can be measured by using 
the reference frame in a distant past. Therefore, a pre- 
diction of expected SN counts as a function of variabil- 
ity flux limits for a given observing time interval, rather 
than as a function of the real SN fiux in a single-epoch 
snapshot, would be useful to quantitatively understand 
the dependence of the event rate on the variability time 
scale. This was not quantitatively presented in earlier 
publications, and here we will present predictions in such 
a form and discuss possible effects of choosing different 
observing intervals. 

The plan of this paper is as follows. In §2 we describe 
basic formulations and calculation methods, and in §3 
we present the results of SN detection rate with some 
discussions on the choice of the observing intervals. In 
§4 we propose the strategy to decipher the CSFH and 
the la delay time, and discuss the feasibility of our strat- 
egy referencing to the plans of future proposed SN sur- 
veys. Discussions and conclusions will be presented in §5. 
Throughout this paper, the standard ACDM universe is 
assumed with following values of the cosmological pa- 
rameters: r^A/ = 0.3, Ha = 0.7, Hq = 70 km s"^ Mpc"^ 
The AB magnitude system is used for the zero-point of 
photometry, unless otherwise stated. 



2. FORMULATIONS AND CALCULATION METHODS 

2.1. Basic Equations 

There are a variety of supernova searches with differ- 
ent search methods and sampling times, but the selection 
process of supernovae can essentially be considered as a 
search of transient objects with a sensitivity magnitude 
limit to fiux variability, miim, between two epochs sepa- 
rated by a time interval T. In reality a more complicated 
selection procedures (e.g., more than two epochs) may be 
taken, but it can be expressed by a combination of multi- 
ple two-epoch selections. The total number of detectable 
SNe, iV(miim, r), in a survey observing an area of flohs 
is calculated from the following equation: 



(iA^i(miim,T, z) 
dz 



where 

rfjVi(TOiim,r, z) 

dz 



n{z) dV{z) 

"obs — J /i(z,TOiim,T) 

1 -\- z dz 



(2) 

Here the subscript i denotes supernova types and the 
summation in Eq. is over all the types of SNe. The 
comoving density of SN rate is described by ri{z), where 
the rate is measured by the cosmic time (or the rest- 
frame time at z). The factor of (1 -I- z)~^ accounts for 
the cosmological time dilation. The comoving volume el- 
ement dV{z)ldz per unit solid angle is in the ordinary 
definition, and fi{z, miim, T) is a function called effective 
visibility time (EVT) or control time, which is defined 
by the observer's time (at z = 0). This function cal- 
culates a time duration in which a SN can be detected, 
and it depends on the spectra and light curves as well 
as the magnitude limit, observing interval, and redshift. 
In a prediction for a snapshot observation without tak- 
ing into account variability, EVT is a duration of a part 
of SN light curve when a supernova is brighter than the 
magnitude limit. On the other hand, in a variability- 
limited survey with an interval of T, the EVT is the 
time duration when the variability fiux is brighter than 
the magnitude limit for variability, i.e., 



i?(mii, 



, ) dt , 



(3) 



where i is a time coordinate for supernova evolution and 
H{x) is the step function [H{x) = 1 and for a; > and 
< 0, respectively], and 



,(z,t,r) = -2.51og{ 
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dz 



(1) 



(4) 

is the magnitude of the variability fiux between two 
epochs of t and t + T. Here m{z, t) is the real magnitude 
of a SN in a given waveband at redshift z and epoch t. 
It should be noted that t is the time for observer and the 
cosmological time dilation must be included in calcula- 
tion of m[z, t). 

2.2. Light Curves and Spectra of SNe 

The apparent SN magnitude m(z,t) can be calculated 
if redshift, light curve, and spectral evolution are known, 
taking into account the K-correction and cosmological 
time dilation by a standard manner. We mostly followed 
the prescription by DF99 for the calculation of supernova 
light curves and spectra, but we collected and used the 
latest data published after 1999 if available, to make the 
most reliable prediction made so far. Here we only de- 
scribe these new data, and other ingredients that are not 
mentioned here are the same with DF99. The lightcurves 
and blackbody temperatures used in our calculation are 
shown in Figure ^ 

The spect r a an d light curves are taken from 
iNueent et al.' ('2002') for type la SNe, while those of SNe 
1999eni (Har auy et al. 2001), 1979c (Branch et al. 1981), 
and 1995G ijPastorello et al.l I2002|) are used for types 
IIP , IIL, and Iln, respectively. The light curves of 
type lb and Ic SNe are similar, and we treat these two 
as a single population, for which the light curve model 
is constr ucted based on the da t a of 19941, 1999ex and 
2002ap dR ichmond et alj Il99(t 'Stritziner er et all l2002t 
iFolev et al. . .2003) . The peak magnitudes differ con- 
siderably from supernova to supernova even within a 
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type, and this is taken into account by integrating the 
peak luminosity function of SNe for all types, assum- 
ing a Gaussian form in magnitude with mean mag- 
nitudes a nd dispersions taken from local values mea- 
sured bv iRichardson et alJ ('2002) (shown in Table [U). 
We use the same relative proportions of the four sub- 
classes of core-collapse SNe with DF99, which are as- 
sumed to be constant with the cosmic time. The sum 
of the relative proportions is not unity because we do 
not include peculiar SN1987A-type supernovae, whose 
peak magnitude is much fainter {Mb ~15) and hence 
they contribute little to detection number of high-z 
SNe. For SNe la, the we ll-known peak-luminosit y and 
lightcurve-shape relation iPhillips 111 993: Phillip s et al.l 



119991: iPerlmutter et all Il997t iStrolger et all l200^ is 

taken into account, while the same shape of lightcurves 
are used for core-collapse SNe for all peak magni- 
tudes. The color /l ightcurve-shape relation of SNe la 
ijPhillips et al.l 1199911 is not included in our calculation. 
We found that this effect would hardly change the results 
presented in our paper. 

We do not take account of the Galactic extinction, 
which can be easily corrected by the known extinction 
map. On the other hand, extinction in host galaxies 
should be dependent on their various properties and 
location of SNe in them, and there should b e a large 
variety llHatano et al.lll99aHTotani & Koba~ shi 1999). 
Since it is difficult to construct a realistic model for this, 
we simply assume a mean reddening of E{B — V) ^ 
0.1 and the standard Milky- Way type extinction curve 
JCardcUi et al. 1989) to make our prediction more real- 
istic than no extinction case. (A detailed discussion on 
this issue will be given in Section 0) 

2.3. the Cosmic Star Formation History 

The CSFH must be specified to predict SN rate as a 
function of rcdshift. Here we use a phenomenological 
CSFH model (Gal- Yam fc Maodl2004 inferred from re- 
cent observations of high-z galaxies: 



\l+z J ^ \ l+z 
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Here a and /3 are the SFR evolutionary indices at high 
(z > 1) and low (z < 1) redshift, respectively. There is a 
consensus that the cosmic SFR increases from z = to 
^ 1-2, whi le the estimated values differ as ^ 1-3 (e.g., 
ILiUv et a ll IT99fit IS uUivan et a ll l2000al: IWilson et al. I 
|2Q02; La nzetta et a l. 2002; Schi minovich et al.ll200,'TH . It 
is even more uncertain whether the SFR density is 
roughly constant to z ^ 5 or declines beyond z 1- 
2 (e.g.. iMadau et al.lll99d: iHuehes et al.lll998(l . In our 
model, a cut off at high redshift z„ = 10 is introduced, 
beyond which there is no star formation, though this cut 
off hardly affects our results. 

The supernova rate density is cal culated by t he usual 
manner as in earlier papers (Madau 'eTal] [19981 DF99). 
The rate of core-collapse SNe (all types except la) is 
simply proportional to the CSFH, but type la rate is 
converted from the CSFH by convolving the delay func- 
tion, i.e., the distribution function of the delay time tjj 
from star formation to the type la events^. We mainly 

^ In some earlier publications, the delay time tjj is defined as the 



use a standard exponential form for the delay func- 
tion, fn(tn) oc eyiv(-tp/T,.) fe.g. IMa,da,u et a1 J 11993 
IGal-Yam fc Ma o?'200l 'Str oker et al J 1200 41) in this pa- 
per. Only in M.ll we calcul ate using the Gaussian type 
function defined iStroker et al.l l)2004D to check the dif- 
ference between them. 

In most earlier papers, the normalization of SN rate 
history is determined by the number of supernova events 
from a unit mass of star formation, but this is highly 
uncertain. To avoid this uncertainty, in this paper we 
always renormalize the ca lculated SN rate history to the 
local SN rate measured bv lCannellaro et all l|1999f) . This 
prescription is also adequate for our purpose, i.e., deci- 
phering CSFH and rja only based on supernova data. 
The renormalization is done separately for type la and 
core-collapse SNe, whose local rate density is 3.0 x 10~^ 
and 7.5 x 10~^ yr~-^ Mpc~^, respectively (|CaDDcllar(J 
12003(1 . Therefore the SN detection number calculated 
in this paper can be considered as a ratio to the local 
SN rates. It should be noted that, with the formulations 
above, the functional form of stellar initial mass function 
(IMF) does not affect the SN rate calculation, provided 
that the IMF is universal and constant with time. 



3. THE SUPERNOVA DETECTION RATE AND OBSERVING 
TIME INTERVALS 

In this section we present some results on the expected 
SN detection rate. To see the dependence on the model 
parameters that we are interested in, we use two CSFH 
models and two values of the la delay time, which are 
largely different but still in the reasonable range. For one 
CSFH model, called the "peak" model, (a,/3) = (-1,3) 
are used, while (a,/3) = (0,1) are used for the other, 
called the "flat" model. The la delay time is set to rja = 
0.3 and 3.0 Gyr, which are referred as "short" and "long" 
models, respectively. 

The left panel of FigI21 shows the total number of de- 
tectable SNe as a function of the magnitude limit for 
variability flux, which is a similar plot to popular galaxy 
counts or star counts as functions of magnitude limit. 
(Therefore we call this plot as "SN counts", hereafter.) 
The SN counts depend on the interval of observation; 
shorter intervals predict less number of detectable SNe. 
The long baseline case (> 360 days) correspond to the 
early studies like DF99, and our results are mostly con- 
sistent with those of DF99. 

The right panel of figure 12 shows the ratio of SN 
counts with different intervals of 30 and 360 days, 
A^(30)/iV(360), for different models of the SN rate his- 
tory. The dependence of this ratio on the SN rate history 
is not large [A{7V(30)/iV(360)} < 0.03]. 

When SNe are selected by flux variability in two epochs 
of observations, SNe can be classified into two categories: 
SNe whose second-epoch flux increases or decreases com- 
pared with the first epoch. We call these as brightening 
and fading SNe. (Therefore it is different from the deriva- 
tive of light curves at a given time.) The ratio of number 
of brightening to fading SNe, Ns/Np, is an observable 

time from the end of the main-sequence phase to the occurrence of 
SNe la. However, such definition is dependent on stellar IMF, and 
here we define it from the time of star formation. For a typical IMF 
(the Salpeter's function with the cut-offs at 0.1 and 100 Mq), the 
mean difference between the two definitions weighted over stellar 
mass is about 0.3 Gyr. 
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quantity that can easily be obtained in any supernova 
surveys, and which is shown in Figure El This ratio is 
determined by the observing interval and shape of SN 
light curves, and potentially it could give some informa- 
tion of SN type mix. When the two images are separated 
enough, both the number of the brightening and fading 
phase should be the same, and the ratio becomes unity. 
For a fixed interval, the ratio becomes smaller at deeper 
magnitude limit, because the faint and long fading phase 
of relatively close supernovae mainly contribute to Np. 
On the other hand, the ratio becomes larger than the 
unity at the brightest limiting magnitude in short inter- 
val observations. In this case only SNe having large vari- 
ability can be detected, and the initial rapid brightening 
phase contributes significantly. Because of the combined 
effects of long decaying phase and cosmological time dila- 
tion, longer intervals are required in deeper observations, 
to get two independent reference frames [Ns/Np = 1]. 
One year interval is enough for this in most cases, but it 
starts to be insufhcient for the deepest surveys possible 
in the near future (Iab ^ 30). 

The right panel of figure 13 shows change of Nb/Np by 
different CSFH models and/or type la delay time. The 
difference is not particularly large and it seems difficult 
to use the observational quantities of A^(30)/A^(360) or 
Nb/Np to discriminate different SN rate histories. How- 
ever, these quantities change significantly with the survey 
depth and time intervals, which can easily be measured 
by observations. Therefore comparing these predictions 
with actual data will provide a useful consistency check, 
such as examination of possible contamination by non- 
SN events like AGNs. 

4. DECIPHERING THE CSFH AND THE la DELAY TIME 

In the framework described above, the three parame- 
ters that we want to constrain by future observations are: 
(1) the low-rcdshift CSFH index, /?, (2) the high-redshift 
CSFH index, a, and (3) the la delay time, Tja- The ques- 
tion that we want to answer is whether we can constrain 
these three parameters by future supernova data, break- 
ing the degeneracy among them. For this purpose we try 
the following strategy. It is expected that, in relatively 
shallow surveys, only supernovae at 2; < 1 are detectable 
and the a parameter does not affect significantly the ob- 
servable quantities, allowing us to constrain mainly /? and 
Tia- (The size of the effect of changing a will quantita- 
tively be checked later.) The survey time interval is set 
to 30 days in all the results shown in this section, and 
our main conclusions do not depend much on the choice 
of the time interval. 

4.1. Determining (3 and ria by SNe ai z < 1 

Then our first approach is to constrain j3 and ria by two 
independent observational quantities at relatively shal- 
low supernova surveys. The quantity that is the most 
easily obtained is the total detected number of SNe in- 
cluding all types at the detection limit of a survey. For 
the second observational quantity, we examine the red- 
shift distribution and the fraction of type la supernovae 
(/la)- Redshift distribution is easier to obtain by pho- 
tometric redshifts of host galaxies without spectroscopy, 
than /la that requires spectroscopy for a firm determi- 
nation. However, we find that there is a considerable 
degeneracy between [3 and Tia in the redshift distribu- 



tion, and the type la fraction has a stronger power to 
constrain /3 and ria- We use a = for the results pre- 
sented in this section. 

4.1.1. Redshift Distribution 

In order to see the power of breaking degeneracy by 
redshift distribution, we calculate the redshift distribu- 
tions and type la fractions using two cases where the val- 
ues of (3 and ria are very different but the total expected 
number of SN detection is the same. Figure^ shows the 
results, where the used parameters are (/3, Tia[Gyr]) — 
(3.0, 4.0) and (2.3, 0.3). It can be seen that there is only 
a small difference in the redshift distribution of SNe of 
all types. The mean redshift differs only by ~ 0.2, and it 
is comparable w ith the accuracy of ph otometric redshifts 
of galaxies fe.g. iFurusawa et alJl200fil) . Even if accurate 
spectroscopic redshifts are available, such a small differ- 
ence of distribution may be hidden by uncertainties in 
theoretical modeling. 

Gal- Yam & Maoz (2004, hereafter GM04) investigated 
how the three parameters (a, /3, ria) can be constrained 
by the redshift distribution of about 1,000 SNe in future 
surveys, and they concluded that the redshift distribu- 
tion gives some useful constraint, which appears some- 
what different from our findings here. In fact, there is no 
discrepancy; we made a similar analysis to GM04, and 
confirmed that our results are mostly similar to those 
in GM04. Here we note that GM04 checked only one 
parameter set, (a,/3, ria) = (—2,4,1), as the input pa- 
rameters for simulations. On the other hand, we found 
the above degeneracy of the redshift distribution by using 
smaller values of (3. This suggests that the degeneracy 
becomes more serious than found by GM04 for smaller 
values of (3. We also note that, even for the input pa- 
rameter set chosen by GM04, the degeneracy between the 
three parameters cannot be completely broken, though 
some parameter space can certainly be ruled out (see Fig. 
7 of GM04). We consider that the redshift distribution 
only is not sufhcient to measure all the three parameters 
by breaking the degeneracy completely. 

4.1.2. The Type la Fraction 

On the other hand, the difference of type la fraction 
is much larger than the redshift distribution in Fig. 01 
indicating that this quantity is a better tool to constrain 
the parameters of interest. For the purpose of seeing this 
clearer, we made contour maps of the expected SN detec- 
tion number and type la fraction, as a function of the two 
parameters of (3 and ria , which is shown in Fig jSl The 
contours of the detectable SN number and la fraction 
are mostly perpendicular to each other, demonstrating 
that they are good sets of observable quantities to break 
the degeneracy. Currently Iab magnitude of 24 is 
the spectroscopic limit for the largest ground-based tele- 
scopes, and it is difficult to measure the SN types beyond 
this magnitude. In this case the expected difference of la 
fraction is not very large (/la ~ 0.5-0.6 for an interest- 
ing range of ria). On the other hand, future space borne 
spectroscopy or extremely large ground-based telescopes 
(30 m class) would allow the type classification at much 
deeper magnitudes. If SN types at Iab ^ 27 are reliably 
determined, we expect that the la fraction differs by a 
factor of about 2 for the same range of ria. 
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In this analysis we have assumed that the effect of a 
on our predictions for shallow surveys is small. If the SN 
la fraction is reliably measured, the parameter /? can be 
estimated without knowing a, since the core-collapse SN 
rate evolution to z ~ 1 simply reflects the SFR evolu- 
tion. However, uncertainty of a affects the estimate of 
Tia. We found that, when the a is varied from —2 to 
1, the estimate of ria by the same observational quanti- 
ties is changed by 1-2 Gyr for the limiting magnitude of 
Iab = 24-27. This degeneracy can be broken by using 
deeper SN counts at magnitude of > 30 (see below). 

4.1.3. On the Form of the la Delay Function 

The exact form of delay function i s unknown, though 
there are s o me theoretical mode ls l|Yungelson fc Liviol 
^M^MM)- IStrolger et all (plfil found that the Gaus- 
sian form of the delay function, foit) oc exp[— (i — 
Tia)2 / {2af^)] with crja ~ (0.2-0.5)Tia, fits better than the 
exponential form used so far in this paper. Therefore we 
repeated the above calculations but changing the delay 
function into the Gaussian form, to estimate the error on 
Tia- It changes by about 1 Gyr, for the li miting magni- 
tude of Iab ~ 27. However, as shown bv'Strolg er et al.l 
12004) ■ the difference of the delay function appears most 
clearly in the redshift distribution of SNe, and hence the 
degeneracy with the delay function form may be removed 
by using the redshift distribution. 

4.2. Determining a by SNe beyond z > I 

It is obvious that we have to observe SNe well beyond 
z ~ 1 to determine a, and detection limit of Iab ~ 27 is 
still not sufficient for this purpose (see Fig. 0)). Future 
projects such as SNAP and JWST can detect supernovae 
down to ~ 30 magnitudes, and their data will be useful 
to constrain a, though it may be difficult to determine 
the SN types at such faintest magnitude. 

We calculate the expected SN counts varying the a 
parameter and the results are shown in Fig. El In this 
calculation we choose a typical parameter set of (/3, ria) 
= (2.0, 1.0). As expected, the change of SN counts by dif- 
ferent a becomes larger for deeper observation. It should 
also be noted that the difference is much larger in the K 
band than / band, since the redshift distribution extends 
farther in redder bands (z ~ 2 for / and z > 5 for K at 
^ 30 mag). In a K-hand observation with miim = 30, 
the fraction of SNe beyond z = 2 is 30-40 % in the case 
of a = 0. 

This is a calculation for fixed values of (/3,ria). By a 
supernova survey down to Iab 27 with secure type 
classification, the parameter (3 can be constrained well 
by using core-collapse SN statistics, while the degeneracy 
between Tia and a remains. It is then possible to break 
the degeneracy by finding a set of (a, rja) satisfying both 
the type la statistics to Iab ~ 27 and the all-type SN 
statistics to ^ 30 magnitudes. 

4.3. Prospects for the Planned Projects 

The current or past superno va surveys f rom ground, 
such as SCP ijPain et alJl200^ and HZT ijTonrv et all 
l200?f>. have reached a depth of about Iab ~ 24. The 
recent GOODS project has achieved a depth of F850LP 
~ 26 thanks to HST observations, but high signal-to- 
noise spectroscopic observation is still limited to Iab ^ 



24 by the largest ground telescopes. As shown above, 
the type la fraction needs to be determined with an ac- 
curacy distinguishing A/ia ~ 0.1 to determine /3 and 
Tia by surveys limited to Iab ^ 24. Requiring sta- 
tistical signal-to-noise of 5, a required number of SNe 
la becomes iV > 1,000, which is much larger than the 
number obtained so far. Howe ver, the ongoing grou nd- 
based proj ects such as SNLS llPritchet et al.ll200^ or 
ESSENCE l|Garnavich et alJ2002HMatheson et alJ2004() 
will greatly increase the number of supernova to this 
magnitude. These surveys could give some useful con- 
straint on P and rja, depending on the accuracy of /la 
and systematic model uncertainties. 

The situation would be improved if we can measure 
the la fraction down to Iab ~ 27, as shown in Fig. |S| 
Spectroscopic type classification to this depth is impossi- 
ble by current ground-based facilities, but classification 
by information obtained by imaging observations, such 
as light-curve shapes, colors, or host galaxy proper ties, 
may still be possible, as tried bv Strol ger et al.l l)2004D for 
the GOODS data. (See also Poznanski et al. 2002.) SNe 
II have very distinctive light-curves (type II-P or II-L), 
and are typically bluer than SNe la at early times. How- 
ever, the reliability of this kind of approach is still rather 
uncertain so far; only about two thirds of the 42 GOODS 
SNe are secure type identifications (Dahlcn et al. 2004). 
A major breakthrough in more distant future will be 
brought by spaceborne spectroscopy, such as SNAP, re- 
liably measuring /la to Iab ^ 27. The SNAP will detect 
and get high signal-to-noise spectra of a few thousands 
of type la supernovae to z ~ 1.7 (Aldcring et al. 2004), 
which corresponds to the peak magnitude of Iab 27 
(see Fig. 0J. Planned extremely large telescopes (~ 30 
m class) will also be useful for type classification of the 
faintest SNe. 

The spacebo rne projects such as SNAP and JWST 
()Panagia'2 003t) will also allow us to break the degeneracy 
between a and ria by SN counts at magnitude of ~ 30. 
As shown above, SNAP will have a sufficient number of 
supernovae but does not probe very high redshift beyond 
z ~ 1 because of the wavelength coverage. The JWST 
could probe higher redshift, but small field of view (14 
square arcminutes) does not allow to detect many super- 
novae at one snapshot. According to our estimate with 
the parameter a = 0.0 , /3 = 3.0 and t = 1.0, JWST wiU 
find about 10 SNe down to K ^ 30 in one snapshot. To 
measure the supernova counts with an accuracy of ~ 10- 
20% to obtain meaningful constraint on a, about 25-100 
supernovae are required. 

5. DISCUSSION AND CONCLUSION 

In this paper we examined the dependence of observa- 
tional quantities obtained by supernova surveys on the 
CSFH and the delay time of SNe la (ria) in detail, and 
tried to identify the best strategy to get meaningful con- 
straints on them by planned future projects. Specifically, 
we chose three parameters of a, P, and ria, and investi- 
gated how well these parameters can be constrained in 
the future based only on supernova survey data, where 
a and P are the indices of the cosmic SFR evolution 
[cx (1 -I- z)^"' at z > 1 or < 1, respectively. 

Relatively shallow surveys {Iab ^ 27) are relevant 
mainly with the parameters of /? and ria. We have shown 
that it is essential to reliably measure the number frac- 
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tion of SNe la (/la) in a survey sensitive to Iab ~ 27, 
in order to measure these parameters. If a ~ 20% dif- 
ference in /la is clearly distinguished for supernovae at 
magnitude of Iab ~ 27, a useful constraints on /? will be 
obtained by core-collapse SN rate evolution. The con- 
straint on Tia is somewhat degenerate with the unknown 
a, but still it can be constrained within an accuracy of 
~ 1-2 Gyr. Some regions of the parameter space of /3 
and Tia are seriously degenerated for the redshift distri- 
bution, and hence we consider that the la fraction is a 
better quantity for this purpose. However, the redshift 
distribution would be useful to constrain the functional 
form of the delay function (e.g., exponential or Gaus- 
sian). 

Type classification down to Iab ~ 27 may be poten- 
tially achieved by ongoing or near future ground-based 
surveys such as SNLS and ESSENCE, but type classifi- 
cation inevitably depends on less reliable methods based 
on imaging information. A more reliable study will be- 
come possible by spaceborne spectroscopy, as proposed 
for SNAP, or future extremely large telescopes on the 
ground. The parameter a can be constrained with break- 
ing the degeneracy with ria, by ultra-deep SN surveys 
down to magnitudes of ~ 30 without SN type classifi- 
cation, which would become possible by future projects 
such as SNAP and JWST. The former has an advantage 
of larger statistics while the latter has an advantage of 
reaching higher redshifts by observing in longer wave- 
length. It should also be noted that, in the above dis- 
cussion, we assumed that the local supernova rates, both 
for type la and core-collapse SNe, are well measured to 
give a good normalization that can be directly compared 
with high redshift observations. Therefore the effort to 
improve the accuracy o f the local supernova rate is also 
highly encouraged fe.g. lAldering et ani2002|) . 

It should be noted that there is a number of simplifi- 
cations in the calculation presented here, e.g., a constant 
IMF, no evolution in supernova light curves or spectra, a 
universal probability of supernova event occurrence per 
unit mass of star formation, and no evolution in the la 
delay time. These uncertainties might induce some bias 
in the CSFH or la delay time inferred from the SN rate 
history. On the other hand, it is possible to get infor- 
mation for these uncertain factors, if we use the CSFH 
determined by independent other methods (such as SFR 
estimates of high-z galaxies) as an external constraint. 

Dust extinction has been incorporated in the same way 
for the SNe la and core-collapse SNe, with a fixed amount 
of extinction and no evolution with z. We adopted this 
simple treatment because it is quite difficult to construct 
a realistic model of dust extinction based on the present 
knowledge about galaxy evolution. However, it is likely 
that extinction has a wide distribution even within a 
galaxy (Hatano, Branch, & Deaton 1998), and mean ex- 
tinction should also evolve with galaxy evolution (Totani 
& Kobayashi 1999). Because the time delay from star for- 
mation to supernova occurrence is different, core-collapse 
SNe co uld be much more heay ily extincted by dust than 
SNe la (jMannucci et alJl2003|) . If the cosmic star forma- 
tion is heavily hidden by dust at z > 2 as inferred from 
submillimetcr observations (e.g., Hughes et al. 1998), 
it might be appropriate to use different values of a for 
type la and core-collapse SNe. These issues should be 
kept in mind when the strategy discussed in this paper 



TABLE 1 

Peak absolute magnitudes (Vega), their dispersion, and 
Relative Proportions of Supernova Types 



Type 


Mb 


a 


f 


la 


-19.13 


0.56=^ 




Ib/c 


-17.71 


1.39 


0.23 


IIP 


-16.67 


1.12 


0.30 


IIL 


-17.70 


0.90 


0.30 


Iln 


-18.82 


0.92 


0.02 


References. 


— Richardson et al.l 


l2()(11) 


for Mb and tr, and 



IDahlen fc FranssonI flSgst) for /. 

''The dispersion of SNe la is uncorrected about the light-curve- 
shape/luminosity relation. 



is actually applied to future data. To see a typical effect 
of changing extinction, we repeated calculations with no 
extinction for SNe la and twice larger extinction for core- 
collapse SNe. The la delay time and /3 are then affected 
by ^1 Gyr and ~0.5, respectively. 

We presented our calculation of expected supernova 
counts as a function of variability flux in a given survey 
time interval, instead of the real flux. We also exam- 
ined quantitatively the change of the detectable number 
of supernovae with time interval, and calculated the ra- 
tio of brightening and fading supernovae between the two 
epochs of the SN selection, which were not quantitatively 
reported in the earlier literature. These are useful for 
planning the observing intervals for the future observa- 
tions. These quantities can easily be measured, and we 
examined a possibility of using these observational infor- 
mation to constrain the interesting parameters of CSFH 
and Tia, but unfortunately they are not particularly use- 
ful. However, comparing these quantities to the theoreti- 
cal predictions would be an important consistency check 
to examine, e.g., a possible contamination of other ob- 
jects such as AGNs. 

We would like to thank an anonymous referee for many 
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the 21st Century COE "Center for Diversity and Univer- 
sality in Physics" and for TT (16740109) from the Min- 
istry of Education, Culture, Sports, Science and Tech- 
nology (MEXT) of Japan. 
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Fig. 1. — Lightcurve and temperature templates for SNe of various types. The left panel is for the lightcurves used in our calculation 
and the right panel is temperature obtained by blackbody-fits to spectra. These are taken from SN 1999emlHamuv et al. 2001), 1979C 
iBranch et al. 1981) and 1995 G |.Pastorello et al. 20 02) for SN types of IIR IIL and Iln. For Ib /Ic. these are constructed based on the data 
of 19941, 1999ex and 2002ap IHiclimond et al.lll99etrStritzinger et alJI2002t H'blev et al. 1120031) . The data for SNe la is from lNugent et alJ 




lAB lAB 



Fig. 2. — Left : Expected number of SNe as a function of limiting magnitude of variability flux, in a field of 1 deg^. Various intervals of 
observation are assumed as indicated in the figure. The peak CSFH model and ti^ = 1.0 Gyr are assumed (see text). Right : The ratio of 
expected SN number with an observing interval of 30 days to that of 360 days. Different line markings correspond to different models of 
CSFH ("peak" or "flat") and the la delay time (0.3 or 3 Gyr), as indicated in the panel. 
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Fig. 3. — Left : The number ratio of brightening to fading SNe between the two epochs of supernova selection. Various intervals of 
observation are assumed as indicated in the figure. The peak CSFH model and -Qa = 1-0 Gyr are assumed. Right : Different line markings 
correspond to different models of CSFH ("peak" or "flat") and the la delay time (0.3 or 3 Gyr), as indicated in the panels. The time 
interval is 30 days. 




Fig. 4. — Redshift distributions of SNe in a survey limited by Iab = 24 and 27, for two models (thick and thin curves) whose total 
expected SN number is almost the same but the CSFH parameters and type la delay time are largely different. The parameters are (/3, 
T]a[Gyr]) = (3.0, 4.0) for thick lines and (2.3,0.3) for thin lines. The solid and dashed lines show the distribution of total number including 
all SN types and that of SNe la, respectively. 
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Fig. 5. — Contour maps of the detectable SN number (thin curves) in a square-degree field and the fraction of SNe la (thick curves), 
for a survey with limiting magnitude of Iab = 24 and 27. The x-axis is the delay time of SNe la (ti^) and the y-axis is the low-redshift 
evolutionary index of the CSFH (/3 in Eq. |HJ. The search interval is assumed to be 30 days. If the interval is 360 days, the SN counts 
become about twice the values indicated in the figure, while the SN la fraction changes very little (at most by -^0.05). 




Fig. 6. — Detectable SN number normalized by the case of a = as a function of limiting magnitude, for all SN types. Various a 
are assumed as indicated in the figure. The used parameters are (/3, T[a[Gyr]) = (2.0, 1.0). Left and right panels are for I and K' band, 
respectively. 
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